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Abstract 


'One  model  for  the  motion  of  n  charged  particles  on  the  x-axis 
leads  to  a  system  of  delay  differential  equations  with  delays 
dependent  on  the  unknown  trajectories.  If  appropriate  past  histories 
of  the  trajectories  are  given,  say  on  [a,0],  then  for  sufficiently 

p 

small  t  0  one  has  a  system  of  n  ordinary  differential  equations 
of  the  form 

y’  *  f(t,y)  with  y(0)  =  yQ  given.  (*) 

The  function  f,  which  involves  the  known  past  histories  of  the 
trajectories,  is  continuous;  so  existence  of  solutions  is  assured. 
However,  f  does  not  satisfy  the  Lipschitz  condition  usually  usecU  i 


re r,  f  does  not  satisfy  the  Lipschitz  condition  usually  usec^  Mz 

{Jit/ 

for  proving  uniqueness.  J 

/  f  f  4 

The  key  new  result  is  that  the  solution  of  (*)  is  unique  <0 

2  ^ 

Lded,  for  some  integer  m  <  n  , 


provided,  for  some  integer  m  <_  n  , 

fjCt,?)  <  1  for  i  =  1,....,  m,  and 

n2  4 

||f(t,£)  -  f(t,n)||  <  K  E  |g  (t-?,)  -  g.Ct-TK)!  +  K  l  |€.  -  nj, 

~  i  =  l  1  1  11  i=m+l  x 

\ 

where  K  >  0  is  constant  and  each  is  a  continuous  function  \ 

of  bounded. variation.  gg 

This  generalized  Lipschitz-type  condition  is  indeed  satisfied  ^ 
in  the  electrodynamics  case.  The  m  components  of  y  which  play 
the  special  role  in  the  above  uniqueness  criterion  are  the  n(n-l)  jfa 
delays  of  the  original  n-body  problem. 

Eventually  one  finds  that  solutions  of  the  original  equations  of 
motion  exist  and  are  unique  as  long  as  no  two  particles  collide. 


\  \ 
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Consider  n  charged  particles  on  the  x-axis;  and  assume  that 
each  is  influenced  only  by  the  retarded  fields  of  the  others,  and 
no  two  particles  collide  during  the  time  considered. 

Let  x^(t)  be  the  position  of  particle  i  at  time  t,  and  let 
c  be  the  speed  of  light.  Then  at  time  t  particle  J  is  influenced 
by  the  electromagnetic  fields  which  were  produced  by  particle  i 
(i  j*  j)  at  some  earlier  instant  t  -  r.^(t).  The  delay  r^j(t) 
represents  the  time  required  for  fields  to  propagate  from  particle 
i  to  particle  j  and  must  satisfy  the  functional  equation 

cr±j(t)  *  |xj(t)  -  x^t  -  r^Ct))!.  (1) 

To  shorten  the  notation,  we  will  sometimes  represent  Xj(t)  and 
r,  ,(t)  by  x,  and  r. ,,  respectively.  Then  Eq.  (1)  becomes 

lj  J 

criJ  =  lxj  "  xi(t  "  rij)l * 

We  shall  also  use  the  symbol 

Vj  =  Vj(t)  =  x’(t)/c 

for  the  velocity  of  particle  j  in  units  of  c.  [If  t  is  an 
endpoint  of  the  interval  of  definition  of  x j ,  then  xj(t)  *s  a 
one-sided  derivative.] 
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The  following  theorem,  proved  in  [1],  lists  some  basic 
properties  of  Eq.  (1). 

i*  Let  xi  and  xj  be  given  differentiable  functions 
with  |  (t )  |  <  c  and  |xj(t)|  <  c  on  [a, 8),  where  ct  <  0  <  8, 
and  with  x^t)  ^  x^(t)  on  [0,0).  Assume  Eq.  (1)  has  a  solution 
when  t  *  0.  (See  Remark  below.) 

Then  Eq.  (1)  has  a  unique  solution  for  all  t  £  [0,0)  and 
(1)  is  equivalent  to 

Crij  =  aij[xj  ”  xi(t  *  ( 1 ' ) 

where 

aij  =  s&n  Cxj(0)  -  x± (0) ] . 

Moreover 

ru  i  ixj  -  *il/2<;- 

Now  assume.  In  addition,  that  Vj  =  x'/c  and  v,  «  x’/c  are 

J  J 

continuous .  Then 

,  vj  "  vi<t  *  r^) 

rij  '  aij  “  v‘(t  -  r^)  for  0  i  t  <  3.  (2) 

Conversely,  if  r±j  is  a  solution  of  (2)  on  [0,8)  and  if  r^CO) 

satisfies  Eq.  (I1)  at  t  =  0,  then  r.,  satisfies  (1>)  on 
[0,6). 

Remarks.  A  sufficient  condition  to  assure  that  Eq.  (l)  has 
a  solution  when  t  -  o  is 

l*i<t)|  <  e  -  |xj(0)  -  X,(0)|/|o|  for  a  <  t  <  0. 

Note  from  (2)  that  d(t  -  r,j)/dt  >0.  So  t  -  r,,(t)  la 
an  increasing  function  of  t.  In  particular,  if  8  <  •,  then 
llmt-*0-  riJ^t^  exists. 


After  these  preliminaries,  we  obtain  the  actual  equation  of 
motion  for  particle  j  as  follows.  Assume  Eq.  (1)  has  a  unique 
solution  for  each  i  t  j,  sum  the  retarded  fields  produced  by 
all  other  particles  (i  /  j),  and  substitute  these  fields 
into  the  Lorentz-Abraham  force  law  for  particle  j.  The  result 
is 


(1  -  v  f?75 


-  I 


Ki.1  “1.1 +  vl(t  -  ri.1) 

p  > 


i*J  i j  ij 


~  v, (t  -  r, ,) 


where  each  K. 


is  a  constant. 


The  "natural"  initial  data  problem  for  the  system  of  functional 


differential  equations  represented  by  (1)  and  (3)  is  as  follows. 


Problem  P.  Let  continuously  differentiable  functions  ...» 

4>n  be  given  on  an  interval  [ot,0]  with  the  properties  that  each 
1 4> -j[  ( t )  |  <  c  on  [a,0]  and  <(>^(0)  ?  <J>j(0)  when  i  ^  J.  A  solution 
of  Problem  P  will  be  a  set  of  continuously  differentiable  functions 
x^,  . . . ,  xn  on  [a, 8)  for  some  8  >  0  such  that 

(a)  x^(t)  =  ^(t)  on  [ct,0]  for  each  i, 

(b)  x^Ct)  /  xj(t)  on  [0,8)  when  i  /  J, 

(c)  |  v± ( t )  |  <  1  on  [0,8)  for  each  i,  where  v.^  =  x|/c, 

(d)  there  exists  a  set  of  functions  {r^J^j  on  [0,8) 
such  that  Eqs.  (1)  and  (3)  are  satisfied  on  [0,8). 


As  yet,  we  have  made  no  assumption  that  the-  "initial 
trajectories",  ...,  $n,  satisfy  any  particular  equations. 

In  the  context  of  the  functional  differential  equations 
presented  here,  it  will  be  quite  natural  to  assume  that  enough 
appropriate  initial  data  is  given  so  that  each  Eq.  (1)  has  a 
solution  at  t  -  0  when  i  ^  J .  (This  assumption 


was  eliminated 


In  the  two-body  problem  by  Travis  [6].) 

The  question  is,  do  we  need  further  smoothness  conditions  on 
in  order  to  establish  the  existence  and  uniqueness  of  a 
solution  of  Problem  P? 

When  this  problem  was  originally  treated  for  the  case  of  two 
particles,  it  was  assumed  that  each  was  Lipschitz  continuous 

[1],  But  later  consideration  of  the  problem  in  three  dimensions 
— involving  equations  of  neutral  type — suggested  that  this  was  an 
unreasonable  hypothesis  [2].  It  seemed  more  realistic  to  require 
each  <♦»’  to  be  just  absolutely  continuous. 

For  the  collinear  motion  of  two  particles,  mere  continuity  of 
each  turned  out  to  be  adequate  [4].  Then,  for  the  case  of  two 

particles  moving  in  three  dimensions  it  was  found  that  absolute 
continuity  is  sufficient  for  local  existence  and  uniqueness  [3]. 

However,  the  proofs  used  in  [4]  and  [3]  do  not  extend  to  the 
case  n  ^  3. 

To  resolve  this  difficulty  the  following  uniqueness  theorem 
for  ordinary  differential  equations  has  recently  been  proved  [5], 

Theorem  2.  The  system  of  p  equations 
y’  *  f(t,y)  with  y ( 0)  =  yQ 

has  a  unique  solution  (locally)  if,  in  some  neighborhood  of  (o,yQ) 

R1+P,  f  is  continuous  and  for  some  integer  m  €  [0,p] 

f^(t»5)  <  1  for  i  *  1,  m 

and  f  satisfies  the  generalized  Lipschitz-type  condition 

m  p 

| |f(t,c)  -  f(t,n) 1 1  <  K  r  I g. ct-e. )  -  g. (t-n.)l  +  K  E  le  -n 

i-i  11  i-m+l  1  i 

where  ||«|l  is  any  norm  in  Rp,  K  >  0  is  a  constant,  and  each 
•  R  *  R  is  continuous  and  is  of  bounded  variation  on  bounded 
sub intervals. 


The  theorem  proved  in  [5]  is  actually  more  general,  but  the 
special  case  above  will  suffice  for  our  present  purpose. 

The  following  (main)  theorem  says  that  the  collinear  n-body 
problem  does  have  a  unique  solution  (until  a  collision  occurs)  even 
assuming  only  absolutely  continuous  initial  velocities — or  slightly 
less. 


Theorem  3*  Let  4*^,  . . . ,  4>n  be  given  functions  on  [ot,0] 
with  each  <J>^  continuous  and  of  bounded  variation.  Assume  that 

(i)  4>j  ( 0 )  ?  4»i(0)  when  j  /  i, 

(ii)  1 4>!(t)  |  <  c  on  [a,0]  for  each  i,  and 

(iii)  Eq.  (1)  has  a  solution  r*j(0)  at  t  =  0  for  each  i  and  j  /  i. 

Then  there  exists  8  >  0  such  that  Problem  P  has  a  unique  solution 

on  [a,8),  and  either  8  «  +00  or  else  for  some  i  and  j  i 

lim  x.(t)  =  lim  x. (t)  —  a  collision. 

t-8-  J  t-*-8- 

Proof .  For  economy  of  notation,  let  r  represent  the  n(n-l)- 
vector-valued  function  with  components  rij • 

Note  that  the  last  paragraph  of  Theorem  1  essentially  lets  us 
consider  Eqs.  (2)  and  (3)  in  place  of  (1)  and  (3).  The  solutions 
of  Eqs.  (2)  and  (3)  which  will  be  of  interest  must  satisfy 


-6- 


Now  define  =  <j>£/c  for  each  i  and  j  /  i.  Then  in  U 

Eqs .  (2)  and  (3)  are  equivalent  to  the  ordinary  differential  equations 


„  .  V1  -  Sl,1(t  -  rlj> 
1J  fflJ  -  %3(t  -  rij) 


_ Vi  7  -  r  K»  al,1  +  g1.1(t  -  -1.1  > 

(1  ‘  'j)3/2  i/j  *±3  °1J  -  Eij(t  -  rij>  ’ 
and  the  initial  conditions  are  r..(0)  *  r?.(0),  v  (0)  =  4>i(0)/c. 

J  J 

Local  existence  and  uniqueness  of  a  solution,  say  for  0  t  <  h, 
follow  from  Theorem  2  with  p  =  n^  and  m  =  n(n-l).  We  have 
thus  generated  the  unique  solution  of  Eqs.  (2)  and  (3)  on  [0,  h). 
Now  let 


'^(t)  on  [a,  0] 

xi(t)  *  ■  t 

1^(0)  +  fQ  cv^s)  ds  on  [0,  h). 

Then  the  resulting  functions  x  . xn  form  the  unique  solution 
of  Problem  P  on  [a,h). 


Next  define 


8  =  sup  {s  >_  0:  a  unique  solution  exists  on  [a,s)>. 

Then  8  >  0  and  a  unique  solution,  x-^,  ...»  xn,  of  Problem  P 
can  be  constructed  on  [a, 8).  Specifically,  for  each  t  <£(0,8), 
define  x^t)  ■  y^t)  where  y1,  . . . ,  yn  is  any  solution  valid 
on  [a,s)  for  some  s  >  t. 

Suppose  (for  contradiction)  that  6  <  00  and  there  exists 
6  >  0  such  that 

|x^(t)  -  Xj(t)|  >_  6  on  [0,8)  for  each  i  and  j  f  i. 


(6) 


-7- 


It  follows  from  Theorem  1  that  each 
rij 1  5/2c • 

Hence  there  exists  a  £,  [0,1)  such  that 

[vi(t  -  )  |  <_  a  on  [0,8)  for  each  i  and  j  /  1. 

Using  these  results,  we  find  from  Eq .  (3)  that 


a  -  v'y 


(1  -  ^ 


is  bounded  on  [0,8)  for  each  J. 


Since  8  <  °°,  this  means  that  for  some  b  £  [0,1) 
f  v j ( t ) |  £  b  on  [0,8)  for  each  j. 

Now,  defining 

x. (8)  =  lim  x. (t)  and  r?.(8)  =  lim  r.,(t), 
t+8-  J  t+8-  J 

we  obtain  functions  x^,  xn  having  the  same  properties  on 

[a, 8]  as  were  originally  assumed  for  (J^,  . .  . ,  <J>n  on  [a,0]. 

Here  r£j(8)  plays  the  role  of  r|j(0).  This  means  one  can 

construct  a  unique  extension  of  the  solution  beyond  t  =  8  (at 

least  for  a  short  interval)  Just  as  we  did  originally  beyond 

t  *  0.  And  that  contradicts  the  definition  of  8. 

So  if  8  <  °°,  it  follows  that  (6)  fails  and  a  collision 

occurs  at  8. 
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